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Abstract 

We study the longitudinal development of hadronic showers in ice for energies 
up to 100 EeV and its implications for radio and optical Cherenkov emission. 
A small fraction of showers induced by primaries of energy above 1 EeV 
is shown to display the characteristic elongation associated to the Landau 
Pomeranchuck Migdal (LPM) effect. The rest look like ordinary showers 
because the highest energy 7r°'s interact instead of decaying in two photons. 
The LPM effect observed in this fraction of the showers is mainly due to n 
and T]' production and decay. We give parameterizations for the total and 
excess charge tracklengths and for the amplitude and angular spread of the 
electric field spectrum in the Cherenkov direction. Implications for neutrino 
detection are briefly addressed. 
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I. INTRODUCTION 



The detection of low energy neutrinos (MeV-GeV) from the Sun and from cosmic ray 
interactions in the atmosphere is now routine in underground detectors which we hope will 
reveal new aspects of the neutrino sector in the Standard Model 0. On the other hand 
the wealth of information obtained with the detection of few MeV neutrinos from supernova 
SN1987A showed that the enormous potential of neutrino astronomy could be reality in the 
near future. Several experiments are being constructed at present to detect neutrinos at 
higher energies, at which the atmospheric neutrino flux is expected to be sufficiently low 
to obtain new information from extraterrestrial sources. They were conceived to measure 
the Cherenkov light emitted by upcoming muons produced in charged current neutrino 
interactions that cross the Earth. The long muon ranges allow an effective volume sufficiently 
large to answer important questions, for example on the intriguing acceleration mechanisms 
in Active Galactic Nuclei (AGN) and Gamma Ray Bursts (GRB) 0. 

It is becoming increasingly apparent that the highest energy neutrinos will play an im- 
portant role. The identification of the highest energy gamma ray sources with Blazars, 
which are interpretated as AGN with jets pointing in our direction, has shifted theorists to 
model acceleration in the jets. These jets are believed to be particle flows with relativistic 
bulk motion and high associated Lorentz factors which naturally make the emission highly 
anisotropic and boost it to very high energies ||. It is also believed that the infrared back- 
ground prevents the TeV gamma ray emission from reaching the Earth except for a few 
nearby sources ||. The possibility is open that particles are accelerated to much higher 
energies getting absorbed in the different backgrounds. Similarly for GRB's there are ac- 
celeration models reaching up to energies much greater than observed. Both phenomena 
are candidates for the intriguing origin of the highest energy cosmic rays |5|||. If this is 
the case, large fluxes of neutrinos can be expected up to close to the highest energies the 
hadrons are accelerated. In any case EeV neutrinos must exist because of the interactions 
of cosmic rays with the microwave photons and their search can answer important issues 
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with respect to these objects or other potential neutrino sources as well as to the origin of 
the highest energy cosmic rays. 

EeV neutrinos have a high enough cross section that prevents them from traversing 
the Earth without interacting. These high energy neutrinos produce extremely high en- 
ergy showers whose Cherenkov light may also allow efficient detection in the underground 
detectors in project or construction. Alternative detection mechanisms become favourable 
for extremely energetic showers. The radio technique relies on the coherent emission of 
Cherenkov light in radio waves by all shower particles when the wavelength of the emission 
is larger than the shower dimensions |7|. Ice is being considered as a good possibility be- 
cause of the long attenuation length for radio waves ||. Similarly the proposed acoustic 
alternative relies on detecting the coherent sound wave which is produced by the shower 
in ocean water ||. The coherence condition makes the signals scale as the square of the 
incident particle energy and this is the reason why these methods become most promising 
at the highest energies. Depending on the actual fluxes the construction of large enough 
detectors sensitive to the expected fluxes may be cost effective provided these alternatives 
are succesfully implemented [[TTJ. Several estimates of the neutrino expected event rates in 
these detectors have been made []TT|-|TB(]. An interesting possibility is the detection of the 
radio pulses from cosmic ray and neutrino interactions on the Moon surface using radiote- 
lescopes on Earth ||14j| . EeV shower behavior in dense media is important for EeV neutrino 



detection, particularly in water (or ice) as chose for most of the proposed detectors. 

Above PeV energies the cross sections for pair production and bremsstrahlung get sup- 
pressed because the characteristic length of the interaction becomes larger than atomic 
spacing and collective effects of the atomic fields have to be considered. This is the Landau 
Pomeranchuck Migdal (LPM) effect ]T3| . Electromagnetic EeV showers have been shown to 
have huge elongations because of the LPM effect |T6|-|T8|| and hadronic showers are expected 



to behave similarly if they contain PeV photons or electrons. 

Here we study the development of EeV hadronic showers in water concentrating on the 
relevant parameters involved in the emission of Cherenkov light and coherent radio pulses. 



We will use many results presented in the study of electromagnetic showers from the same 
perspective as described in ref [|Tj|. As it was the case for electromagnetic showers, the full 
simulation of Cherenkov light and in particular coherent effects between different particles 
is not practical because of the high energies involved. We choose to approximate shower 
development by a combination of Monte Carlo methods and parameterizations of showers of 
lower energies. The main characteristics of the radio pulse emission from these showers at 
the Cherenkov angle can be deduced from the longitudinal shower development so we will 
use one dimensional approximations which will only be valid for emission in this direction. 
Although more detailed calculations will be necessary, many of the most important aspects 
of the possibilities of detecting hadronic showers as induced by neutrino interactions can be 
answered with the approximations made here. Many of the qualitative conclusions obtained 
in this work can be easily extended to other media. 



II. HADRONIC SHOWERS. 

Shower development has been studied in air up to 10 21 eV because all high energy cosmic 
ray experiments look for them. At lower energies (< 1 TeV) they have also been studied in 
dense media for accelerator experiments. The study of purely electromagnetic showers and 



LPM implications, has been done in water for energies up to 10 EeV | lq| . Such showers are 
produced at the lepton vertex in charged current Deep Inelastic Scattering (DIS) electron 
neutrino interactions. In these interactions the average energy transfer to the lepton is 
predicted to be large, about 75% 

Hadronic showers produced in both charged and neutral current DIS neutrino interactions 
are started by hadrons in the fragmentation of the nuclear debris. At high energies the 
multiplicity is expected to be about half of that in e + e~ collisions for a center of mass 
energy squared equal to minus the square of the 4- momentum transfer (Q 2 = —q^q^) in 
the DIS process |[20|| . Practically all models for neutrino production involve the decay of 
pions which naturally predict more muon flavor neutrinos from naive channel counting. The 
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ratio of muon to electron neutrino production may be actually larger than expected in some 
environments because of muon synchrotron losses as has been recently stressed ||21|| . All 
these factors may well compensate the fact that in DIS interactions only one fourth of the 
neutrino energy is typically transferred to the hadron vertex. Multiple hadrons are also 
produced in 68% of the M^-boson decay in the resonant electron antineutrino scattering 
with atomic electrons. Hadronic showers will also play an important role in tau neutrino 
detection through the proposed double bang events fl22| . 

For EeV neutrinos the fragmented particles, mostly pions, will mantain the direction 
of the primary neutrino in the lab frame because their average transverse momentum is 
expected to be in the few hundred MeV range. These particles continue to interact with 
nucleons and they still have sufficient energy so that after three or four generations the 
cummulative angular deviations are still very small. The hadronic shower can be thought as 
having a hard penetrating central core containing mostly pions which feeds electromagnetic 
subshowers fundamentally through tt° decay in two photons. Because the medium is dense, 
charged pions are expected to interact before decaying. Assuming energy equipartition 
between all flavor pions, a fraction of 1/3 would go into electromagnetic subshowers every 
time there is an interaction. If all interactions occur simultaneously every interaction length 
(to make things simple) it is easy to find that after three generations 80% of the energy is 
in electromagnetic subshowers. The shower will have a markedly electromagnetic character 
because of the high energy involved and the medium density which prevents decays of most 
of the charged pions. The lateral distribution should be related to a sum of distributions 
of electromagnetic showers started at different depths. Near the core the showers are thus 
younger, dominated by small electromagnetic showers initiated deeper, while away from the 
shower axis they are due to early electromagnetic subshowers of higher energies. 

Continuing with this simple model, we can consider that after a neutrino interaction 
with pion multiplicity N, the pion energy is roughly the energy transferred to the hadron 
yE divided by N. The LPM effect is known to increase the electron and photon interaction 
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lengths above some threshold energy ~ Elpm (2 PeV for ice) f] []17 , |16 , ^3[ . However the 



longitudinal profile of electromagnetic showers in ice around shower maximum starts to be 



affected typically at primary energies above E = 20 PeV [|18[]. For a 1 EeV energy transfer 
to the hadron a multiplicity of order 15 can be expected in the v interaction and the photons 
from the 7r° decays would roughly have 30 PeV. Moreover most of the electromagnetic energy 
in the shower is produced after the first generation. As a result little deviations from typical 
shower development are expected below 1 EeV because of the LPM effect. 

As the primary energy rises so will the energy of the first generation 7r°'s. However 
they do not necesarily produce high energy photons. Neutral pions have a decay length 
cr = 2.5 10~ 6 cm and an interaction depth of about 130 gcm~ 2 so that in the competition 
between decay and interaction the latter dominates for ir° energies above about 6.7 PeV. 
As a result the showers can be expected to show LPM effects in a mitigated form. Even 
some 100 EeV showers are actually shown below to display small LPM elongations. EeV 
hadronic showers can look more like rescaled versions of lower energy showers while EeV 
electromagnetic showers are greatly distorted and elongated because of the LPM effect. 

III. SIMULATIONS. 

A complete simulation of the radio pulse should monitor particles down to the MeV 
range requiring unreasonably large times for EeV energies [24j. To investigate EeV hadronic 



shower characteristics we have developed a fast hybrid Monte Carlo which simulates 1 di- 
mensional showers down to some crossover energy, at which the subshower produced is taken 
from a tested parameterization. For electromagnetic subshowers we use the Greisen parame- 



terization which was shown in |24| to be valid up to energies of 100 TeV. We have calculated 
our own parameterizations for hadronic subshowers inititated by protons, pions and kaons. 
The code has been designed to calculate quantities such as the projected tracklength (onto 



''Corresponds to the definition in Ref. [17| 
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shower axis) and the weighted projected tracklength (excess tracklength due to electrons). 
These quantities are respectively known to be the relevant parameters for optical and radio 
Cherenkov emission EM. We have used parameterizations of these tracklengths for electro- 



magnetic showers using the results obtained in ||24|| . Intermediate results will be presented 
elsewhere. 

The shower simulation part is based on UNICAS [2l| which has been adapted to run in a 
homogeneous water medium with density 0.924 gem -3 . The code has also been modifyed to 
take into account interactions of 7r°'s and other shortlived resonances as well as corrections to 
bremsstrahlung and pair production accounting for the Landau- Pomeranchuk-Migdal effect. 
For the lower energy showers we have performed full simulations. Above 100 TeV we use the 
hadronic code SIBYLL together with the parameterizations. SIBYLL includes "minijet 



events" which are responsible for the cross section increase with energy |27]] and soft collisons 
based on the Dual Parton Model (DPM) pSf. 

We have made extensive independent checks of the codes used. For example full 1-D 
simulations (for energies below 100 TeV) are well within the uncertainties associated to 
the different hadronic models compared with CORSIKA in [29]. We have also checked the 
1-D version of UNICAS in ice comparing photon showers with those obtained in the 3-D 



electromagnetic code developed specifically for calculating radio pulses in ice [p4|| . We have 
checked the hybrid approach both comparing it to full simulations and ensuring that all 
hybrid codes are stable under changes of the energy crossover. 



IV. RESULTS. 

We have simulated showers initiated by protons, and most of the results presented here 
correspond to such showers. Little difference is observed between proton and pion initiated 
showers of the same primary energy. We have explored the energy range from 10 GeV to 100 
EeV simulating large numbers of proton showers (over 20,000 showers above 1 EeV) with 
the described codes. The lower energy showers are necessary to make parameterizations 



7 



which allow the simulation of the highest energy ones. 

We firstly show the results of the fraction of energy going into electromagnetic showers 
in Fig. 1. Here the fraction is seen to rise with energy according to the simplified model 
above, until the effect of neutral pion interactions becomes important. In the total track- 
length calculation, relevant for the Cherenkov light output, the electromagnetic contribution 
accurately scales with the electromagnetic energy in the shower as shown in Fig. 1. Other 
charged particles contribute some percentage which decreases with energy (15 % at most) as 
can be expected from the electromagnetic energy fraction in Fig. 1. The coherent radio pulse 
emission scales with the difference between positive and negative charge tracklengths, the 
weighted tracklength, to which particles other than electrons and positrons contribute neg- 
ligibly. We also display "shower length" in Fig. 1, defined in as the length along which 



shower size exceeds 70% of its maximum. The angular width of the Cherenkov radiopulse 
emitted at the Cherenkov angle is approximately inversely proportional to this parameter. 

Hadronic and electromagnetic showers have well known differences which have been stud- 
ied at high energies only in air, mainly hadronic showers have larger attenuation lengthf], 
reach deeper from first interaction point, have harder lateral distributions and carry a sig- 
nificant fraction of their energy in muons and neutrinos. As a result they have slightly less 
particles at maximum than electromagnetic showers of the same primary energy. These dif- 
ferences can be seen in Fig. 2 where proton showers are compared to pure electromagnetic 
showers of energies 10 TeV, 1 and 100 PeV. As long as the LPM effect is not significant, 
the showers "scale" with primary energy in much the same way as electromagnetic showers, 
what incidentally, simplifies their parameterizations. 

Below 1 EeV the longitudinal profile of hadronic showers has this scaling behavior, char- 
acterisitc of showers that are not affected by the LPM effect. This is not surprising, the 



^ Attenuation length is defined as the length over which the energy drops by a factor e in the 
region well after shower maximum and it is proportional to the particle interaction length. 
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average energy of the pions is well below the primary interaction energy because of the large 
multiplicities involved |3(J and only electromagnetic showers of energy above ~ 20 PeV show 
"scaling deviations" [IB]. In table I we give the mean energy of the 7r°'s as a function of 
primary energy as obtained in the hadronic interaction generator. In Fig. 3 we compare 
the energy distributions of all the first generation 7r°'s in proton-ice interactions at different 
energies to those that actually decay into photons. Only about 10% of the 7r°'s are expected 
to produce photons of energy above 20 PeV because most of them interact. Showers are not 
elongated despite being induced by primaries exceeding Elpm by several orders of magni- 
tude. In Fig. 1 we also display shower length neglecting neutral pion interactions to further 
illustrate this point. 

In this respect the role of resonances with shorter lifetimes than the ir° can be more 
important if they have decay channels involving electrons or photons. If these resonances 
are produced in early interactions, high energy photons (or electrons) produce subshowers 
with characteristic LPM elongations. The extent of this effect is a probabilistic issue and 
depends on detailed aspects of the hadronic interaction model. At the highest energies, 
experimental evidence is non-existant and our results have to be taken with caution. They 
indicate that these aspects are important for shower development and it may be that, if these 
showers are ever observed, they would provide experimental information on such channels. 
In particular we find that the r\ and the rj' contribute most and this is just because they are 
very short lived, they are produced relatively often and they have considerable branchings 
into photons, ~ 40% (~ 65%) for r] (r/fl). 

We find that a proportion of showers above 1 EeV have deep tails characteristic of LPM 
showers (see Fig. 4). We also note that the energy content of these tails (proportional to 
shower development integral) is typically of order 1%. All of these can be quantitatively 
addressed with the aid of Fig. 5, representing the "leading photon" distribution as obtained 



^Including the rf — ► rj it + tt decay chain. 
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in SYBILL, that is the differential probability that the highest energy photon in a shower 
has a fraction energy (x = E^/E ) in the range (x, x+dx). The median is at about 2% 
and not very sensitive to primary energy, hence about 50% of the hadronic showers having 
energy above 50 E t are expected to have a photon of energy > E t . If E t is chosen say above 
100 PeV, the electromagnetic subshower generated by the leading photon will have a long 
LPM tail. The range of leading photon energy fractions is between 0.1% and 30% in 99.9% 
of the showers and one can easily deduce that it is practically impossible to have a shower 
without an LPM tail if the primary energy is above 100 EeV and that hadronic showers of 
energy below 250 PeV will hardly have LPM effects. 



V. CHERENKOV RADIATION FROM HADRONIC SHOWERS. 

Cherenkov light is emitted as charged particles travel through a medium at a speed 
greater than that of light. While the electric field amplitude is proportional to particle 
tracklength, in the optical region the emission is incoherent and the power output scales 
with the total tracklength |H| . For wavelengths exceeding the shower dimensions Cherenkov 
light is emitted coherently by all shower particles, opposite charges contribute with opposite 



phases and the electric field amplitude is proportional to the weighted tracklength p2 



The excess electrons in electromagnetic showers are caused by low energy interactions with 



atomic electrons [2Jj3^j. The matter antimatter symmetry in the cross sections is broken by 
the absence of target positrons. The dominant process is Compton scattering with atomic 
electrons, but electron positron annihilation, Moller and Bhabha scattering also contribute 



to an overall 21% excess when expresssed in terms of the total tracklength p4 |. Many of 
the results obtained in the context of radio emission from electromagnetic showers fl8| , |2"4}| 
can be applied to our hadronic shower simulation results. 

There is a smooth transition from the incoherent to the coherent Cherenkov regions 



which has been described in |[18||2^| . It occurs in ice at ~1 GHz for radiation emitted in the 



Cherenkov angle direction and at lower frequencies away from this direction. The spectrum 



10 



rises linearly with frequency until interference effects from different shower regions start 
to take place. This happens at a maximum frequency v max . For pulses observed at the 
Cherenkov angle the interference is governed by the lateral distribution, the narrower the 
lateral spread the bigger the maximum frequency. We expect it not to be too different from 
electromagnetic showers, and we take v max from them. As a first approximation this is 
justified because most energy in hadronic showers is electromagnetic, and consists of much 
lower energy subshowers. 

The electric field amplitude in the Cherenkov direction can be approximated quite ac- 
curately as a 1-D Fourier transform of the longitudinal charge excess distribution. The half 
width of the central peak is inversely proportional to shower length as defined above and 
shown in Fig. 1. This has been used to describe Cherenkov emission from EeV electromag- 
netic showers in ice ||18|| . Since the highest output signal is concentrated in the Cherenkov 



peak, its size and angular spread are the most relevant aspects of such radiation, particularly 
for establishing the detection capabilities of future detectors. These properties are calculable 
in this approach and can be used for detector design while more detailed simulations are 
implemented. 

The change in shower length with energy is very small, ~ 30% between 100 TeV and 
10 EeV as shown in Fig. 1. As a result the angular distribution of the Fourier transform of the 
radiopulse is not going to become narrower for higher energies contrary to electromagnetic 
showers |18|] . For those showers that do display LPM effects the narrowing can be observed in 



a rather convoluted diffraction pattern, see Fig. 6; the Cherenkov radiopulse being coherent 
contains all the information on the excess charge distribution. The lateral distributions of 
the excess charge also affects the full diffraction pattern behavior in other directions but this 
is not addressed in this work. 

From our shower results it is straigthforward to obtain a parameterization for the total 
tracklength t from which the total Cherenkov light in the optical band can be obtained: 

/(e) = -1.27 KT 2 - 4.76 KT 2 (e + 3) - 2.07 10~ 3 (e + 3) 2 + 0.52 v / eT3 (1) 
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e = 10910 1^ (2) 

t = 6.25 /(e) m (3) 

The weighted projected tracklength is ~ 0.21 t and the normalization of the radiopulse 
electric field spectrum at the Cherenkov peak E(u, R, 9 C ) scales with it: 

R\E{u>,R,9 c )\ = LI x 10- 7 f(z) ^ L i + .4(^ V MHZ_1 (4) 

where 1/ is the observation frequency and z/ ~ 500 MHz. 

Lastly the angular spread of the electric field is inversely proportional to shower length: 

-ln2 



E(u, R, 9) = E(u, R, 9 C ) e L A9 J (5) 



' 1° ^ (2.07 - 0.33 e + 7.5 10" 2 e 2 ) for 1 TeV < E < 100 TeV 

1° f (1.74 - 1.21 10- 2 e) for 100 TeV < E < 100 PeV (6) 
k 1° f (4.23 - 0.785 e + 5.5 10" 2 e 2 ) for 100 PeV < E < 10 EeV 



VI. DISCUSSION. 

We have shown that EeV hadronic showers produced in neutrino interactions are very 
different from electromagnetic showers. The discussion above applies to showers initiated 
by protons or pions of a given primary energy. A 100 EeV neutrino transferres 25% of its 
energy to the hadronic debris with typical multiplicity 17 and only the leading particle keeps 
a significant fraction 1 — Koi the primary neutrino energy, where K is the mean inelasticiy. 
As a result practically all of the electromagnetic subshowers will have energies far removed 
from the neutrino energy. We can expect the hadronic showers induced by neutrinos up to 
10 3 EeV to have a quite ordinary longitudinal development curves around shower maximum. 
They will however differ considerably from the electromagnetic showers produced in charged 
current electron neutrino interactions which display a typical fluctuating nature, reaching a 



great deal further and having many less particles at shower maximum |T8 
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It should be in principle possible to distinguish them comparing the shower profiles 
if they occur inside optical Cherenkov detectors such as those now in construction. For 
radiodetection the differences are converted to differences in the angular width around the 
Cherenkov direction. This effect favours the detection of hadronic induced pulses because 
the antenna spacing would not have to be reduced. It is not clear that at EeV energies 
the electron neutrino contributions for the electromagnetic pulses dominate. The answer 
depends crucially on the relative electron to muon neutrino fluxes, the fraction of energy 
transferred to the nucleons in DIS interactions and the relative spacing between antennas. 
With an array of antennas sufficiently dense it should be possible to recognize the much 
broader pulses for EeV hadronic showers. This is very interesting since the relative rates of 
each type can be related to the ratio of electron to muon neutrinos. This could be used to 
test neutrino mixing parameters models in regions of parameter space unaccesible to other 



experiments and similar to those suggested in Ref. |33 
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FIGURE CAPTIONS 



Figure 1: Left: Electromagnetic energy fraction as a function of primary energy (solid) 
and ratio of the weighted projected tracklength and shower energy in units of km TeV -1 . 
Right: Shower length (in radiation lengths and defined in text) in a proton shower as a func- 
tion of primary energy. The solid (dashed) line represents shower length when competition 
between interaction and decay for 7r°'s is (is not) performed. 

Figure 2: Shower development of a proton (solid) and electromagnetic (dashed) showers 
of primary energies 10 TeV, 1 PeV and 100 PeV. 

Figure 3: Energy distribution of 7r°'s in 10 EeV, 1 EeV and 100 PeV proton-ice inter- 
actions showing the average number of 7r°'s with energy greater than xE n (10 4 events). The 
three overlapping upper curves correspond to all the 7r°'s produced while the three lower 
curves count only those 7r°'s that decay. 

Figure 4: Shower development of individual showers of energies 10 EeV (solid his- 
tograms) and 1 EeV (dashed lines). For 1 EeV we also show the average over 100 showers 
(dotted line). Labels refer to the individual showers taken to obtain Fig. 6. 

Figure 5: Leading photon distributions for 10 EeV and 100 EeV hadronic showers. The 
area under this graph, as shown, is proportional to the probability. 

Figure 6: Radiopulse distribution around the Cherenkov direction for the two showers 
labelled in Fig. 4, illustrating the effect of the LPM tail on the radiopulse angular distribu- 
tion. 



TABLES 



E (TeV) 


10 4 


10 5 


10 6 


10 7 


< E n o > (TeV) 


40 


300 


2 10 3 


1.7 10 4 



Table I: Average neutral pion energy in proton proton collisions for different projectile 
energies in the lab frame as produced by SIBYLL hadronic generator. 



17 



x/N dN/dx 




Mod(E(o;)) (arbitrary units) 




